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The influence of large amplitude pairing fluctuations is investigated in the framework of beyond mean field symmetry 
conserving configuration mixing calculations. In the numerical application the finite range density dependent Gogny 
force is used. We investigate the nucleus 54 Cr with particle number and angular momentum projected wave functions 
considering the axial quadrupole deformation and the pairing gap degree of freedom as generator coordinates. We 
find that the effects of the pairing fluctuations increase with the excitation energy and the angular momentum. The 
self-consistency in the determination of the basis states plays an important role. 
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It is well known that pairing correlations play a pre- 
dominant role in most nuclear phenomena [TJ [2] . In spite 
of that nuclear superfluidity is not completely understood 
and it is still a very active research field [3j 2J [5j [6] . Nowa- 
days pairing correlations are routinely included in all mean 
field nuclear structure calculations either in the Hartree- 
Fock (HF) plus BCS approach or in the Hartree-Fock- 
Bogoliubov (HFB) approach. However while the HFB 
(BCS) approach is very successful to describe infinite sys- 
tems, it was soon recognized [2] that it is a deficient theory 
in the case of atomic nuclei where only a few pairs of nu- 
cleons (weak pairing regime) participate in the superfluid 
phase. Furthermore in a finite system a sharp number of 
particles is needed. 

In spite of this fact the modern beyond mean field the- 
ories (BMFT) that have been developed in recent years 
[7] do not treat pairing on the same footing as the nuclear 
deformation. From those theories the so called Symme- 
try Conserving Configuration Mixing (SCCM) approaches 
have played a major role. Their two basic ingredients 
are: first the restoration of the symmetries broken at the 
mean field level, i.e., the angular momentum projection 
-and sometimes particle number projection- and the sec- 
ond one the configuration mixing. The latter one is done 
in the Generator Coordinate Method (GCM) which allows 
to deal easily with collective degrees of freedom. The dif- 
ficulty to solve the corresponding variational equations, 
however, increases considerably with the number of co- 
ordinates and until now only the deformation parameter 
f3 has been considered in axially symmetric calculations 
|9j [TDJ [11] and (/?, 7) and K mixing in triaxial ones 
[I2j [I3j [14] . It is the purpose of this letter to consider 



explicitly for the first time the collective pairing degree of 
freedom with a finite range force in the framework of the 
SCCJVQ approach with axially symmetric shapes. 

Some arguments indicating the relevance of considering 
the pairing fluctuations are: First, the monopole (pairing) 
and the quadrupole (deformation) are the most relevant 
degrees of freedom, it seems therefore reasonable to con- 
sider both of them on an equal footing. Second, the con- 
strained variational principle of Ritz used to determine the 
intrinsic basis states is very effective in determining the 
wave function (w.f.) of the ground state with the given 
quantum numbers and constraints. Ground states of the 
SCCM calculations benefit of this fact, however, SCCM ex- 
cited states with the same or different quantum numbers 
are not favored by it and depend more strongly on the basis 
size (number of generator coordinates). Consequently, in 
restricted self-consistent calculations a stretched spectrum 
is forecasted which will be squeezed by an appropriate in- 
crease of the basis size, for example by allowing pairing 
fluctuations. Lastly, the pairing vibrations on their own 
are very interesting and a simultaneous study on shape 
and pairing fluctuations will allow us to disentangle if there 
exist genuine pairing vibrations or alternatively to which 
degree they are damped. 

In the gauge space associated with pairing the HFB 
w.f. has two collective degrees of freedom, the pairing gap 
A, which measures the amount of pairing correlations, i.e., 
the "deformation" [15] in the associated gauge space, and 
the angle cp which indicates the orientation of the HFB 
state in this space. The HFB minimization determines 
the w.f. and thereby A while the gauge angle cp does not 
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play any role at the mean field level. The degree of free- 
dom associated to ip has been exploited in the past [16] : 
linear combinations of w.f.'s with different orientation in 
the gauge space provide a number conserving wave func- 
tion. Pairing vibrations -associated with w.f.'s with differ- 
ent pairing gaps- around the average gap parameter Aq of 
the energy minimum, on the other hand, have attracted 
little attention. As a matter of fact they have been con- 
sidered only either with very schematic interactions in the 
framework of the collective Hamiltonian [T71 HE] , hi micro- 
scopic model calculations [19] [20], in reduced configuration 
space [21 or in earlier BMFT approaches [221123]. 

To quantify the pairing content of a w.f. with a fi- 
nite range interaction, like the Gogny force, that provides 
state dependent gaps is not trivial. A quantity that sup- 
plies a measure of the pairing correlations and is easy to 
handle is the mean square deviation of the particle num- 
ber (A7V) 2 . This quantity is zero in the absence of pair- 
ing correlations and is large for strongly correlated sys- 
tems. Furthermore, since for a schematic pairing interac- 
tion <(AiV) 2 ) = 4 £ fc>0 u\v\ = A 2 £ fc>0 ^ with E k the 

k 

quasiparticle energy, A ex ((A7V) 2 ) 1 / 2 and (A7V) 2 provides 
an indication of the pairing content of the wave function. 
In the following; we will denote S = (^(ATV) 2 ^) 1 / 2 and 
use it as coordinate to generate wave functions with dif- 
ferent pairing correlations. Of course one could constrain 
separately the mean square deviation for protons and for 
neutrons, this is more complicated and work in this di- 
rection is on progress. In this first exploratory letter we 
constrain only the total mean square deviation. 

In order to implement pairing fluctuations together 
with axially symmetric quadrupole fluctuations we pro- 
ceed in the following way: First, we generate intrinsic 
HFB wave functions \(/>(q,6)) with given quadrupole de- 
formation q and "pairing deformation" 5 by solving the 
variational equation 

5E'[<p(q,5)] = 0, (1) 

with 

E' = i ^^ L - XMQM ~ \sm^N?W'\ (2) 

and the Lagrange multipliers X q and being determined 
by the constraints 

(<t>\Q2o\<t>) = q, m&NfW 2 =8, (3) 

Q20 is the quadrupole mass operator. If in Eq.J^J |<I>) = \<j>) 
we are solving the plain HFB equation^] A Particle Num- 
ber Projection (PNP) and/or Angular Momentum Projec- 
tion (AMP) out of this w.f. would be a Projection After 
Variation (PAV). However if |$) = P N \</>), being P N the 



2 In this case we have to add to Eq. [2]a term —XjyN to keep the 
particle number right on the average, with Ajv fixed by the constraint 

(c/>\N\c/>) = N. 



particle number projector, the determination of |</>) is done 
in the so called Variation After (Particle Number) Projec- 
tion (PN-VAP) approach. This method provides a much 
better description of the pairing correlations in the intrin- 
sic w.f. [24] although is more involved. There are only 
few implementations of the PN-VAP approach using ei- 
ther separable forces [25], small configuration spaces [26] 
or the most recent ones with the Gogny [27] and Skyrme 
functionals [28, 29 j. Finally, as in the HFB case, an an- 
gular momentum projection can be performed afterwards. 
The variational equations are solved using the conjugate 
gradient method [30]. Once we have generated the basis 
states we can proceed with the configuration mixing calcu- 
lation. This is carried out within the generator coordinate 
method taking linear combinations of the wave functions 
obtained in the first step after performing projections on 
the required symmetries 

= J f'^q, S) P \<P(q, 5))dqdS. (4) 

The variational principle applied to the weights f I,(J (q,S) 
gives the generalized eigenvalue problem, the Hill- Wheeler 
(HW) equation: 

J (H 1 ^, q'S') - E^AT^qS, q'S')) (q' '^dq 1 'dS' = 0, (5) 

with H 1 and Af 1 the Hamiltonian and norm overlaps, re- 
spectively, see [9 for further details. The symbol a = 
1,2,..., numbers the ground and excited states with angu- 
lar momentum J, E I,a is the energy of the corresponding 
state. If in Eq. [4] the projector P = P z ' P N P 1 , with P 1 a 
short-hand notation for the angular momentum projection 
operator, we are performing a PNP and AMP. In previous 
calculations [9] [10] one can find configuration mixing cal- 
culations with only AMP, in this case P = P 1 . We shall 
also discuss below such kind of calculations. 

We use in our numerical application the finite range 
density dependent Gogny force with the DIS parametriza- 
tion [31], which is well known for its successful predictions. 
In the configuration space we take eight oscillator shells; 
for the q coordinate we take an interval —220 fm 2 up to 
+400 fm 2 with Aq = 20 fm 2 and for the pairing coordinate 
5 from to 4.5 with a 0.5 step size. The chosen 5 interval 
covers a pairing energy range from 0.0 up to ~ 50 MeV, 
to compare with typical values of a few MeV's in the nu- 
cleus we shall analyze. In both cases the intervals have 
been chosen in such a way that an energy convergence is 
obtained and that the tails of the collective wave functions 
go to zero inside the interval border. 

As an example we will apply our theory to the nu- 
cleus 54 Cr. Our first intention is to investigate a numer- 
ical case to quantify the relevance of the pairing fluctu- 
ations and not to make an exhaustive comparison with 
the experiment. We are aware that relevant degrees of 
freedoms like triaxiality, or the use of a larger configura- 
tion space to describe very excited states, are necessary 
for such a goal. A comparison with the experiment for the 
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Figure 1: (Color online) Left: Potential energy curves as a function of 
the quadrupole moment for the nucleus 54 Cr in various approaches. 
The energies of (a) have been shifted by 20 MeV. Right: The wave 
functions of the lowest 0+ collective states (0^~ continuous line, 0^" 
dashed line and 0^!" dotted line ) in the corresponding approaches. 
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Figure 2: Contour plots of the square root of the absolute value of 
the pairing energy in MeV in the (q, 5) plane. 



2 + states in this region without pairing fluctuations can 
be found in ref. [TT]. In order to disentangle the effects of 
particle number projection we will present three different 
methods, namely, HFB+AMP, HFB+PNAMP and PN- 
VAP+PNAMP to determine the w.f.'s P\<j)(q, S)) of Eq.jiJ 
see Table [U 

In the first one (HFB+AMP) particle number projec- 
tion is completely ignored, i.e., the intrinsic w.f., 
is determined in the HFB approach and afterwards AMP 
is performed. In the second one (HFB+PNAMP) | (/>(<?, 5)) 
is determined in the HFB approach but then PNP and 
AMP are performed. In the last one (PN-VAP+PNAMP) 
\4>(q,5)) is determined in the PN-VAP method and after- 
wards the PNP and AMP are performed. This w.f.'s can 
be used by their own (to calculate potential energy surfaces 
for example) or as basis states of the SCCM calculations, 



Name of the 
Approach 


Intrinsic w.f. \<p) fixed 
by minimization of: 


p 


HFB+AMP 




P 1 


HFB+PNAMP 


(4>\H\<f>)/(<f>\<f>) 


pIpN 


PN-VAP+PNAMP 


{<t>\HP N \4>)/{4>\P N \4>) 


pIpN 



Table 1: Approaches used to calculate the w.f. P\cj){q, 8)) of Eq. [4] 
In the name of each approach the letters before the plus sign indicate 
the way in which the intrinsic w.f. \(/>(q, 5)) is determined. The letters 
after the plus sign indicate the symmetries we are considering in the 
projection P. 



i.e., to solve Eq.|5j we will present examples of both cases. 
See also Tabl^l] Furthermore in order to investigate the 
effect of the pairing fluctuations we shall present results 
in one dimension (ID), i.e., only with the quadrupole mo- 
ment as coordinate generator and in two dimensions (2D) 
with the pairing energy and the quadrupole moment as 
generator coordinates. 

Prior to the SCCM calculations we present in Fig.[TJa) 
an example of potential energy curves as a function of 
the quadrupole moment (ID) calculated within the three 
approximations described above for the nucleus 54 Cr. All 
the three calculations present two minima, one prolate and 
one oblate, being the HFB energy (dotted line) the flattest 
one. The HFB+PNP wave functions (dash-dotted line) 
produces an additional energy lowering with the excep- 
tion of some values around q ~ 120 fm 2 where the pairing 
correlations are zero. This is not the case in the PN-VAP 
approach where a larger energy gain is obtained for all q 
values. In Fig.[TJb) we display the AMP curves for / = h 
where we can observe the effects of the angular momentum 
projection (additional energy gain and slight displacement 
of the minima to higher deformations). Now we proceed 
to the SCCM approaches: In Fig. [TJc)-(e) we show the 
wave functions of the ground and first two excited states 
for I = Oh obtained by solving Eq.|5]for the different cases. 
These collective wave functions represent the distribution 
of probability of having the states with a given intrinsic 
deformation^) We discuss the SCCM results for the PN- 
VAP+PNAMP case as an example. We observe in Fig 1(e) 
that the structure of the states is strongly correlated to the 
shape of the potential energy curve. The 0+ state presents 
two peaks at q = 120 fm 2 (the largest) and q = — 80 fm 2 , 
the Oj state has also two peaks but this time the oblate 
one is the largest one. Finally, the O3" state presents a nar- 
row peak structure at q = 240 fm 2 corresponding to the 
shoulder of the potential energy curve at this deformation. 

We now discuss the inclusion of the 6 degree of free- 
dom. In Fig. [2] we present contour lines of the square root 
of the pairing energy in the (g, S) plane for the wave func- 
tions of the PN-VAP+PNAMP approach for I = Oh. We 
obtain almost straight lines indicating that, as in the BCS 



3 See [2] for the relationship of the weights G used to plot the w.f. 
and the weights / of Eq. ^] 
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Figure 3: (Color online) Energy contour plots in the plane (8, q) in different approaches for 54 Cr, see main body for further explanations. The 
dashed lines represent contours from to 3 MeV in lMeV steps. The continuous lines represent contours from 4 to 10 MeV in 2 MeV steps. 
The energy origin has been chosen independently for each panel and the energy of the minimum has been set zero. The energy scale on the 
right is in MeV. 



case with schematic interactions, the pairing energy is pro- 
portional to (AN 2 ) justifying the use of the constraint on 
(AN 2 ) 1 / 2 to generate wave functions with different pairing 
content. As mentioned above we can see that the pairing 
energy interval covered by our calculations ranges from 
zero to about 50 MeV and by looking at this figure one 
can have an idea of the behavior of the pairing energy at 
each point of the plane (q, S) in this case. 

As with the one dimensional case prior to the SCCM 
results we first present in Fig. [3] contour lines of the poten- 
tial energy of 54 Cr as a function of the constrained param- 
eters (g, S) in different approaches. The bullets represent 
the S values of the self-consistent solution (HFB or PN- 
VAP) extracted from the one dimensional (g-constrained) 
approach and are displayed as a discussion guide. Since 
all HFB based approaches do have the same intrinsic w.f. 
all of them have the same bullets pattern. The same ap- 
ply to all PN-VAP based approaches. We first discuss the 
HFB based approaches. In the panel (1,1), here we use 
matrix element notation (file, column), we display the two 
dimensional potential energy in the plain HFB approach, 
i.e., the solution of Eqs. |1|2| with |<£) = \<j>). Since in 
this case the bullets represent the self consistent solutions 
the line of the bullets is, as expected, orthogonal to the 
equipotential contour lines. As a matter of fact following 



the bullets we can reproduce the corresponding curve in 
Fig. [TJa). The gross behavior is that for |g|-values larger 
than 50 fm 2 the equipotentials are parallel (perpendicu- 
lar) to the J-axis for 5- values smaller (larger) than a given 
^max(^)- That means, for a given q we can modify the 
pairing content of the wave function up to a given 5 max (q) 
without much energy cost, i.e., the Hartree-Fock average 
potential is able to readjust itself to compensate for the 
pairing energy gain, producing very soft potentials in the 
S direction. However, for S values larger than £) max , rather 
suddenly, this compensation does not take place anymore 
and it costs a large amount of energy to increase the pair- 
ing correlations of the system, producing very stiff poten- 
tials. For |g|-values smaller than 50 fm 2 , i.e., close to the 
spherical shape, there is a valley along the value 5 ~ 2.5, 
i.e., the energy increases if one tries to increase or decrease 
the pairing correlations. This behavior has probably to do 
with the fact that we are close to the spherical shape and 
dealing with very pure configurations: four protons in the 
f?/2 and two neutrons in the sub-shell. This constel- 
lation corresponds to a given pairing energy, and since the 
other sub-shells are far away in energy it takes energy to 
increase or decrease the pairing correlations. 

We can get an idea about the effect of the PNP just 
taking the HFB w.f.'s of the previous calculations and eval- 
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HFB 


HFB+PNP 


PN-VAP 


Basic approach 


ID 


-470.096 


-470.708 


-473.066 




2D 


-470.096 


-471.620 


-473.066 


AMP (1 = 0+) 


ID 


-472.490 


-472.903 


-475.860 




2D 


-472.491 


-474.535 


-475.860 


AMP (1 = 2+) 


ID 


-472.017 


-472.348 


-475.023 




2D 


-472.017 


-473.473 


-475.025 


SCCM (0+) 


ID 


-473.522 


-474.985 


-476.636 




2D 


-474.137 


-475.809 


-476.865 


SCCM (0+) 


ID 


-471.282 


-470.285 


-471.566 




2D 


-473.182 


-470.720 


-472.232 


SCCM (2+) 


ID 


-472.541 


-473.639 


-475.407 




2D 


-473.096 


-474.260 


-475.598 



Table 2: Energy minima, in MeV, in different approaches in ID and 
2D. In the first row the values obtained in the basic approaches (HFB, 
etc) . In the second row the energy obtained performing the AMP to 
/ = + of the corresponding column, i.e., the number in the "HFB" 
column corresponds to the HFB+AMP approach, the one in the 
column "HFB+PNP" to the HFB+PNAMP y the one in'PV-VAP" 
to the PN-VAP+PNAMP. In the third row the same as in the second 
one but for 7 = 2+. In the fourth row the energy of the solution of 
the Hill- Wheeler equation jS]), (SCCM) for the 0+ state taking as 
basis states the corresponding to each column as explained for the 
second row. In the fifth and sixth rows the same as in fourth row 
but for the 0^ and the 2+ states, respectively. 

uating the PNP energy in each (g, S) point. This is the 
HFB+PNP approach which we display in panel (1,2) of 
Fig. [3J As compared with the plain HFB results we find 
that the minima are somewhat deeper and slightly shifted 
to spherical shapes. One also observes a small shift of the 
equipotentials to larger 5 values. In particular the abso- 
lute minimum has been shifted to a value of S ~ 2, in 
contrast to the plain HFB solution. We can also compare 
with the one dimensional plot of Fig. [TJa). We see now 
in the two-dimensional plot that the kink that appears 
around q ~ 120 fm 2 will not be there if we shall follow 
the one-dimensional path perpendicular to the equipoten- 
tials. In panel (1,3) we present the results of the PN-VAP 
approach. Here we observe an enhancement of the effects 
mentioned with respect to the HFB+PNP plot, namely 
a shift of the equipotentials to higher S values and more 
pronounced minima. It is also worthwhile to observe that 
the ID solution in this approach represented by the bullets 
(see also Fig. [TJa)) goes perpendicular to the equipoten- 
tials in the 2D surface. 

We now discuss the effect of the AMP (for / = Oh) on 
top of the discussed approximations. That means we take 
now the intrinsic w.f. determined in the HFB and the PN- 
VAP and calculate the AMP and/or PNP energy in the 
plane (q,S). In panel (2,1) we display the HFB+AMP re- 
sults. As compared with the HFB approach ( see panel 
(1,1) of Fig. [3J, we find a softening of the equipotentials 
for large q values and the reinforcement of both minima. 
Interestingly the AMP on top of the HFB does not change 



the position of the minima in the (<?, S) plane. In panel 

(2.2) we present the HFB+PNAMP results, as compared 
with the HFB+PNP approach (panel (1,2)), we observe 
a deepening and a shift of the minima to larger q- values. 
However, we do not see a shift in the 5 direction. In panel 

(2.3) the PN-VAP+PNAMP results are shown. They dif- 
fer from the PN-VAP ones in a drift of the equipotentials 
to larger deformations and a deepening of the minima. 
Concerning the 5 coordinate the positions of the minima 
are the same as without AMP. A comparison between the 
PN-VAP+PNAMP approach and the HFB+PNAMP one 
reveals that the minima of the latter one are much softer 
towards small 5 values than the former ones. Finally in 
the third row of Fig. [3] we present the results of the PN- 
VAP+PNAMP approach for higher angular momenta. As 
the angular momentum increases we observe a weakening 
as well as a shift towards sphericity of the oblate mini- 
mum and a drift of the equipotentials towards smaller 5 
values. This effect shows clearly that with growing angular 
momentum smaller pairing energies are preferred. 

Let us now consider the results of the SCCM calcula- 
tions, for which we have to solve the HW equation (Eq.[5|. 
We again analyze the three cases we have just discussed. 
To evaluate the impact of the pairing fluctuations on the 
different observable we also consider the solutions of the 
HW equation in ID with one, (<?), and in 2D with two, 
(#,£), generator coordinates. We have calculated the four 
lowest states for each angular momentum and ordered the 
levels according to their quadrupole moment. Concern- 
ing the ground states of the three approaches the middle 
panels of Fig. [3] are a graphical representation of the di- 
agonal matrix elements of Eq. [5] We will use these panels 
to guide the analysis of the final results. However, since 
non-diagonal matrix elements do no follow general sim- 
ple rules it is not always easy to understand fine details 
of the spectra. The magnitude of the non-diagonal ele- 
ments and thereby the energy gain after the solution of 
the HW equation depends on the approach. In general, 
because of self-consistency and the quality of the approach 
before the HW diagonalization we expect the smallest en- 
ergy gain in the PN-VAP+PNAMP case, a larger one in 
the HFB+PNAMP and the largest one in the HFB+AMP 
approach. In Table [2] we provide the energy values at 
the minima in different cases. In the 2D case the en- 
ergy lowering of the 0^ states in the SCCM approach 
(4th row) with respect to the HFB+AMP (2nd row) is 1.0 
MeV (PN-VAP+PNAMP), 1.27 MeV (HFB+PNAMP) 
and 1.65 MeV (HFB+AMP). 

In Fig. [4] we display the spectra in the different ap- 
proaches. We start with the SCCM results for the PN- 
VAP+PNAMP case. When we compare the one dimen- 
sional (ID) with the two dimensional (2D) spectra we ob- 
serve that the general behavior is to compress the other- 
wise too much stretched ID spectrum. The general ten- 
dency is that the lowering of the states increases with the 
excitation energy and with the angular momentum. The 
reason for this behavior, as mentioned above, is the vari- 
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Figure 4: (Color online) Spectra of 54 Cr in different approaches. In the bottom of each spectrum the approach used to generate the 
Hill-Wheeler basis is indicated. Dashed lines (quantum numbers at the right hand side) correspond to ID calculations with the quadrupole 
moment as generator coordinate. Continuous lines (quantum numbers at the left hand side) correspond to 2D calculations with q and 5 as 
generator coordinates. 



ational principle. To solve the HW equation in ID the 
intrinsic w.f.'s of the bullets displayed in panel (1,3) of 
Fig. [3] are needed, and for the 2D the ones of the bul- 
lets plus all the others. Because of self-consistency the ID 
path of the bullets goes along minima and saddle points. 
One should also notice that in the frame of the cranking 
model the condition (J x ) = + 1) for the case I = ft 

coincides with the variational principle of Eq. [I] This is 
the reason why the bullets "almost" provide the right path 
along the PN-VAP+PNAMP potential energy surface for 
/ = h in panel (2,3) of Fig. [3j One should notice, for 
example, that for 7 = 6, see panel (3,3) this is not the 
case anymore. However since the Yrast band is rather ro- 
tational the energy gain of the I = 6± state in the 2D 
case as compared with the ID is not as large as for 6^, 63" 
or 6 4 states. Based on the comments above we expect 
the deviation of the 2D and ID results to be such that 
the smallest energy gain takes place for the ground state. 
Furthermore, larger deviations occur for non- Yrast states 
hih etc, in such a way that the larger the angular mo- 
mentum and the energy, the larger the deviation. We now 
consider the SCCM spectra for the HFB+PNAMP case. 
In this case the ID results resemble the PN-VAP+PNAMP 
ones to some extend. The 2D results however behave in a 
different way as compared to the PN-VAP+PNAMP. The 



main discrepancy is the fact that most of the 2D states 
lie below the corresponding ID ones and in particular the 
ground state band is more stretched in 2D than in ID. The 
reason for this is the lack of self-consistency clearly seen in 
panels (1,2) and (2,2). Here we observe that the bullets do 
not proceed through the minima and one cannot take for 
granted that the ground state is the one with the lowest 
energy gain, see Table |2j Furthermore since in the ID cal- 
culations the pairing energy is small at the minimum the 
moment of inertia is large producing a more compressed 
spectrum than the 2D in which the minima have stronger 
pairing correlations. 

Our last case is the SCCM results for the HFB+AMP 
approach. The first observation is the general compres- 
sion of the ID and the 2D calculations. If we compare 
the three ID calculations we find that HFB+AMP are the 
more squeezed ones. The reason can be easily seen com- 
paring the panels (2,1), (2,2) and (2,3) of Fig.||and Fig.[l] 
There we can see that the minima in the HFB+AMP solu- 
tions, i.e., the place where the collective w.f.'s concentrate 
do have much less pairing correlations than in the minima 
of the other approaches, i.e., larger moments of inertia. 
In particular, we can also compare in Table [3] the pairing 
energy of the 0^ states in the ID approaches. When we 
include the pairing fluctuation, since we are dealing with 
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Figure 5: (Color online) Contour plots of the SCCM collective wave functions in the (£, q) plane of the three lowest 1 = states in the different 
approaches. In the insets of each plot the approach used to generate the Hill-Wheeler basis is indicated. 



self-consistent calculations, though in a smaller proportion 
than in the case of the PN-VAP+PNAMP, we expect a 
compression of the spectrum. The compression, however, 
is huge as compared with that case and the HFB+PNAMP 
one. Since the only difference between the HFB+AMP 
and the HFB+PNAMP approaches is the PNP we con- 
clude that this ingredient is very relevant. The fact that 
the Yrast band is more compressed in the ID than in the 
2D at variance with the PN-VAP+PNAMP approach is 
obviously related to the magnitude of the non-diagonal 
matrix elements as compared with the diagonal ones. In 
Table [3] we can observe larger pairing correlations for the 
ground state in the 2D than in the ID calculations. 

We can have a clearer interpretation of the results look- 
ing at the collective wave functions. In Fig. [5] we display 
the collective wave functions of the 0^,0^ and O3" states 
in the three theoretical approaches. The SCCM w.f's of 
the PN-VAP+PNAMP case displayed in the upper pan- 
els, are easy to understand. Looking at the panel (2,3) of 
Fig. [3] we find two minima and a softening of the contour 
lines at q « 240 fm 2 . At these positions we find maxima 
in our wave functions. The 0^ state has a strong peak 
on the prolate side which is rather soft in the 5 direction 
and a smaller one on the oblate side. The 0^ state has 
a peak on the oblate side, with a two hump structure in 
the 5 degree of freedom, and a weaker one on the pro- 
late side. The O3" state has a big peak at a rather high q 



value and two smaller ones. As a matter of fact a compar- 
ison with the ID results in the right panel of Fig. [I] shows 
that the collective 2D w.f.'s peaks at q values close to the 
corresponding ID ones. The main characteristic of the 
collective PN-VAP+PNAMP w.f.'s is that they are very 
soft in the S direction. Concerning the HFB+PNAMP 
w.f.'s, see second row of Fig. [5j they have the general ap- 
pearance as the PN-VAP+PNAMP though they are not 
that soft in the S direction. The O3" state is presenting 
more mixing on the prolate side than the corresponding 
PN-VAP+PNAMP one. The big difference, as expected, 
appears in the behavior of the HFB+AMP collective w.f.'s. 
We find the strength more concentrated than in the other 
two cases and we find much more mixing. In particular 
the prolate and oblate peaks of the 0+ and 0^ states in 
the PNP approaches do mix in the HFB+AMP leading 
to a concentration of probability at small prolate shapes. 
The same applies to the O3" state, we find an unnatural 
mixing of the different minima giving rise to a strange ac- 
cumulation of the strength providing the collapse of the 
spectrum as we have found in Fig. [4j These results could 
be an indication that the configuration mixing calculations 
with w.f.'s belonging to weak and strong pairing regimes 
do not perform properly if we treat the particle number 
symmetry at first order of the Kamlah approach. 

As mentioned above variational approaches are spe- 
cially suited for the description of the ground state. We 
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2+ 


2+ 

z 3 


Ez En 


Ez En 


Ez En 


Q spec 


Q ] spec 


Q ] spec 


PN-VAP+PNAMP-1D 


-4.8 -5.4 


-4.8 -4.9 


-4.4 -5.5 


-29.670 


16.845 


-49.563 


PN-VAP+PNAMP-2D 


-4.9 -5.6 


-4.1 -3.9 


-4.8 -6.9 


-29.556 


15.299 


-42.722 


HFB+PNAMP-1D 


-2.2 -2.2 


-2.0 -2.6 


-2.6 -3.3 


-26.997 


14.631 


-34.864 


HFB+PNAMP-2D 


-3.6 -5.0 


-2.5 -2.8 


-2.9 -5.7 


-27.950 


13.824 


-40.985 


HFB+AMP-1D 


-1.3 -2.0 


-1.8 -1.7 


-0.2 -3.2 


-30.184 


-6.766 


-27.401 


HFB+AMP-2D 


-1.7 -3.0 


-3.3 -4.1 


-2.9 -5.5 


-31.224 


-20.381 


-26.925 



Table 3: Proton and neutron pairing energies, in MeV, of the first 0+ states of the SCCM calculations as well as the spectroscopic quadrupole 
moments, in e 2 fm 2 , of the first 2+ states. In each row the approach used to generate the Hill-Wheeler basis is indicated. 



can study the convergence of the ground state energy with 
the improving approaches. As we can see in Table [2j 
the energy in the basic approach, i.e., the plain HFB, is 
—470.096 MeV. If we now consider self-consistently the 
particle number symmetry, i.e., the PN-VAP approach, 
we obtain a lowering in the energy of 2.97 MeV. The pro- 
jection of the PN-VAP wave function to zero angular mo- 
mentum provides an additional decrease of 2.79 MeV. The 
mixing of the quadrupole degree of freedom gives a further 
lowering of 0.78 MeV and, finally, the additional considera- 
tion of the pairing fluctuations contributes with 0.23 MeV. 
These results indicate a good convergence of the ground 
state energy. 

One could also rise the question about pairing vibra- 
tions. Looking to the collective w.f. of the PNP ap- 
proaches we do not find clear evidence of genuine pairing 
vibrations. What one observes is a strong softening of the 
peaks of the w.f.'s with the pairing degree of freedom. In 
Table [3] we present the pairing energies of the three lowest 
+ states in the three approaches and in the ID and 2D 
cases. In the PN-VAP+PNAMP case the largest values 
are obtained and no big difference between the ID and 2D 
cases is observed. Smaller absolute values are found in the 
HFB+PNAMP approach, a large increase in the neutron 
pairing energies is observed for the 0+ and O3" states. In 
the HFB+AMP approach one observes that the pairing 
energies in the O3" are much larger in the 2D case than in 
the ID case. In a boson theory [15] a pairing vibration is 
interpreted as a two boson state and one could be tempted 
to conclude that this state is of such character. This nu- 
cleus is probably not a good candidate to look for pairing 
vibrations and further work must be done in other regions 
of the nuclide table. 

The spectroscopic quadrupole moments of the three 
lowest 2 + states are also shown in Table [H In the SCCM 
approaches with PNP the 2+ and 2^ states are prolate 
and the 2^ oblate, both approaches providing comparable 
results though the PN-VAP+PNAMP ones are somewhat 
larger. The SCCM approach with HFB+AMP w.f.'s pre- 
dicts prolate deformation for the three states and larger 
deviation from the PN-VAP+PNAMP one. In general for 
this nucleus we do not find big changes between the ID 
and 2D quadrupole moments. 



In conclusion, for the first time we have considered si- 
multaneously pairing and quadrupole fluctuations in the 
framework of the symmetry conserving configuration mix- 
ing approach with effective forces. In the studied nucleus, 
54 Cr, we find a large effect of the pairing fluctuations on 
the energies of the excited states. We find that particle 
number projection is very relevant to obtain physical re- 
sults. We also obtain a strong dependence on the way the 
intrinsic basis is generated, the variation after projection 
approach providing the most consistent way. The spectra 
obtained without particle number projection provide very 
strong mixing causing an unreasonable compression of the 
spectrum when pairing fluctuations are included. In this 
nucleus no clear evidence of genuine pairing vibrations is 
found but a strong softening of the quadrupole modes with 
the pairing degrees of freedom. Neither we find a strong 
influence of the pairing fluctuations on the quadrupole mo- 
ments. 
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1. Appendix 

As a guide to the reader we put together of all the 
acronyms in the order of appearance in the text. 

HF Hartree Fock 

BCS Bardeen Cooper Schrieffer 

HFB Hartree Fock Bogoliubov 

BMFT Beyond Mean Field Theory 

SCCM Symmetry Conserving Configuration Mixing 
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GCM Generator Coordinate Method 
PNP Particle Number Projection 
AMP Angular Momentum Projection 
PAV Projection After Variation 
VAP Variation After Projection 

PN-VAP Particle Number Variation After Projection 



References 

[1] A. Bohr, B. R. Mottelson and D. Pines Phys. Rev. 110 (1958) 
936. 

[2] P. Ring and P. Schuck, The Nuclear Many Body Problem 

(Springer- Ver lag, Berlin, 1980). 
[3] T. Lesinski, T. Duguet, K. Bennaceur, J. Meyer, Eur. Phys. J. 

A 121 (2009) 

[4] G. Potel et al, Phys. Rev. Lett. 105 (2010)172502 

[5] N. Pillet, N. Sandulescu, P. Schuck, and J.-F. Berger, Phys. 

Rev. C 81 (2010) 034307 
[6] A. Gezerlis, G. F. Bertsch, and Y. L. Luo Phys. Rev. Lett. 106 

(2011)252502 

[7] M. Bender, P.-H. Heenen and P.-G. Reinhard, Rev. Mod. Phys. 
75 (2003) 121 

[8] M. Bender, P. Bonche, and P.-H. Heenen, Phys. Rev. C 74 
(2006)024312 

[9] R. Rodriguez-Guzman, J. L. Egido, L. M. Robledo, Nucl. Phys. 
A 709 (2002)201 

[10] T. Niksic, D. Vretenar, and P. Ring, Phys. Rev. C 74 

(2006) 064309 

[11] T.R. Rodriguez and J.L. Egido, Phys. Rev. Lett. 99 

(2007) 062501 

[12] M. Bender and P.-H. Heenen, Phys. Rev. C 78 (2008)024309 
[13] J. M. Yao, J. Meng, P. Ring, and D. Vretenar, Phys. Rev. C 81 
(2010)044311 

[14] T.R. Rodriguez and J.L. Egido, Phys. Rev. C 81 (2010)064323 
[15] R.A. Broglia, O. Hansen and C. Riedel Advances in Nuclear 

Physics, Vol. 6 (M. Baranger and E. Vogt, Eds.), Plenum Press, 

New York (1973 

[16] K. Dietrich, H. J. Mang and J. H. Pradal, Phys. Rev. 135 (1964) 
B22. 

[17] D. Bes, R. Broglia, R. Perazzo and K. Kumar, Nucl. Phys. A143 
(1967) 1 

[18] A. Gozdz, K. Pomorski, M. Brack et al Nucl. Phys. A 
442(1985)50 

[19] M. A. Fernandez and J.L. Egido, Phys. Rev. B68 (2003)184505 
[20] M. A. Fernandez and J.L. Egido, Eur. Phys. J. B 48 (2005) 
305317 

[21] A. Faessler et al Nucl. Phys. A 217(1973) 420 
[22] J. Meyer et al., Nucl. Phys A533(1991)307 
[23] P.-H. Heenen et al, Eur. Phys. J. All(2001)393. 
[24] M. Anguiano, J. L. Egido, and L. M. Robledo, Phys. Lett. B545 
(2002) 62 

[25] J. L. Egido and P. Ring, Nucl. Phys. A 383 (1982) 189; Nucl. 

Phys. A 388 (1982) 19. 
[26] K. W. Schmid and F. Gruemmer, Rep. Prog. Phys. 50 (1987) 

731. 

[27] M. Anguiano, J.L. Egido and L. M. Robledo, Nucl. Phys. A 

696 (2001)467 
[28] M.V. Stoitsov et al Phys. Rev. C 76 (2007)014308 
[29] G. Hupin, D. Lacroix, M. Bender, [arXiv: 1105 .4084 and private 

communication. 

[30] J.L. Egido, J. Lessing, V. Martin, L.M. Robledo, Nucl. Phys. A 
594 (1995) 70. 

[31] J. F. Berger, M. Girod, D. Gogny, Nucl. Phys. A 428, 23c (1984) 



9 



